SOLITON INTERACTION WITH SLOWLY VARYING POTENTIALS 



JUSTIN HOLMER AND MACIEJ ZWORSKI 



Abstract. We study the Gross-Pitaevskii equation with a slowly varying smooth poten- 
tial, V(x) = W(hx). We show that up to time log(l/ h)/h and errors of size h 2 in H l , 
the solution is a soliton evolving according to the classical dynamics of a natural effective 
Hamiltonian, (£ 2 + sech 2 * V(x))/2. This provides an improvement (h — > h 2 ) compared to 
previous works, and is strikingly confirmed by numerical simulations sec FigJTJ 

1. Introduction 

The Gross-Pitaevskii equation is the nonlinear Schrddinger equation with an external 
potential: 

id t u + \d 2 u — V[x)u + u\u\ 2 = 
u(x, 0) = e lv ° x sech(x - a ) , 



In [9] we investigated the case of V(x) = h 2 W(x) where < h < 1 and W G i/ _1 (R; R). 
Because of our earlier work on high velocity scattering [7], [8], which is potential specific, the 
paper presented the case of W(x) = 5 (x), but as was pointed out there the method applies 
verbatim to the more general case, W G H~ x . Motivated by the approach of [H] we now 
revisit the case of slowly varying potentials, V = W(hx), and show that using the effective 
Hamiltonian approach we can describe the evolution of the soliton with errors of size h? . In 
particular, in this setting, we improve the results of Frohlich-Gustafson-Jonsson-Sigal [1]. 

Theorem 1. Suppose that in (Ol) . V{x) = W(hx), where W G C 3 (M; E). Let 5 G (0, 1/2). 
Then on the time interval 

8log(l/h) 

we have 

(1.2) \\u(t,m) - e ^W e i7(i) sech (._ a ( t ))|| H1(R) < Ch 2 ~\ 0<h<h , 

where a, v, 7, solve the following system of equations 



1 

a = v , v = — 



-sech 2 * V'(a) 



:i-3) , 



2 



7 = - + — — sech 2 * V(a) + (x sech 2 x tanh x) * V(a) , 

with initial data (a ,v ,0). The constants C and h depend on \\W^\\loo , < k < 3, and 
\v \ only. 
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Figure 1 . Comparison of the dynamics of the center of motion of the soliton 
for the Gross-Pitaevskii equation with a slowly varying potential, 

iut = — -u xx — \u\ 2 u — sech 2 (hx)u , h—1/5, ft, = 1/4, 

and initial condition in (11.11) with vq = 0, = —3. The dashed red curve 
shows the solution to Newton's equations used in [1] and [I], the blue curve 
shows the center of the approximate soliton u, and the black dashed curve is 
given by the equations of motion of the effective Hamiltonian 

^ (v 2 + sech 2 (/i«) * sech 2 (a)) . 

The improvement of the approximation given by the effective Hamiltonian 
is remarkable even in the case of h — 1/4 in which we already see radiative 
dissipation in the first cycle. 
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As in [9] the proof of our theorem follows the long tradition of the study of stability of 
solitons which started with the work of M.I. Weinstein [H]. The interaction of solitons 
with external potentials was studied in the stationary semiclassical setting by Floer and A. 
Weinstein [3] and Oh pj2|, and the first dynamical result belongs to Bronski and Jerrard 
PP, see also [2], [II]- The semiclassical regime is equivalent to considering slowly varying 
potentials, and the dynamics in that case was studied in [I] , [2] , [B] , [13] (see also numerous 
references given there). 

The results of [I] in the special case of (11. ip give 
(1.4) |K*,«) -e^ w e i7(t) sech(«-a(t))|| H i (M) < Ch 1 ^ , 0<t< 5\og{l/h)/h, 
where 

a = v + 0(h 2 ), v = -V'{a) + 0(h 2 ), 
(1 ' 5) 7 = 1/2 + v 2 /2 - V(a) + 0{h 2 ) , 

with initial data (do, vq, 0). In other words, the motion of the soliton is approximately given 
by Newton's equations, a — v, v — —V'(a). However it is not clear if, as in (jl.3p . the errors 
0{h 2 ) in (11. 5p can be removed without affecting (11.41) - see $61 Since (I1.3P imply equations 
(II. 5p . Theorem [H shows that we can replace h x ~ 5 by h 2 ~ 5 in (jl.4p . keeping (II. 5p . 

When V(x) = h 2 W(x), and W G [9], Newton's equations are clearly not 

applicable. To describe evolution up to time log(l/ 'h)/h we introduced a natural effec- 
tive Hamiltonian. A numerical experiment shown in [§1 Fig. 2] and reproduced here in 
Figfj] suggested that the effective Hamiltonian approach gives a dramatic improvement for 
slowly varying potentials. Theorem [U (and a more precise Theorem [2] below) quantify that 
improvement by changing the errors from h in (11.41) to h 2 in (11.21) . 

To describe the natural effective Hamiltonian we recall that the Gross-Pitaevski equation 
(II. ip is the equation for the Hamiltonian flow of 

(1-6) H v {u) d ^ \ J {\d x u\ 2 - \u\ A )dx + l -J V \u\ 2 , 

with respect to the symplectic form on i/ 1 (M, C) (considered as a real Hilbert space): 

(1.7) u(u, v) = Im J uv , u,v E if 1 (M, C) . 

When V = 0, r] = sech is a minimizer of H Q with the prescribed L 2 norm (||r/||| 2 = 2): 

(1.8) d£ v = 0, £{u)^ H Q {u) + \\\u\\ 2 L2 . 

Here dE n is the differential of 8 : H 1 — > R, see §2.11 
The flow of H is tangent to the manifold of solitons, 

M = {e i7 e™ (x - a) /isech(/i(x-a)), a,w, 7 GM, fi G R+} , 
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which of course corresponds to the fact that the solution of (11.11) with V = and u (x, 0) = 
e i70+ ™ o(x - ao) /isech(/i(x - a )), is 

U(x 7 t) = g(t) ■ 7] = e H+i«o(a : -ao)+i(/, 2 -^)*/2 ^ sec h( At (x - Oq - V t)) , 

(1-9) df 

g (t) = (a + uo*, «o, 7o + (At 2 + ^ 2 )/2, fJ.) . 

The symplectic form (11.71) restricted to M is 

(1.10) (jJ\m = ^dv A da + vdfi A da + c?7 A d[i , 

see §2.31 The evolution of the parameters (a, v, 7, /i) in the solution (11.91) follows the 
Hamilton flow of 

H \ M - — --, 

with respect to the symplectic form uj\m- 

The system of equations (11.31) is obtained using the following basic idea: if a Hamilton 
flow of H, with initial condition on a symplectic submanifold, M, stays close to M, then 
the flow is close to the Hamilton flow of H\m- In our case M is the manifold of solitons 
and H is given by fl 1 . 6 [) 

2 3 1 

(1.11) H V \ M (a, v,j,n) = f^--^- + -fi 2 (V * sech 2 )^) . 

2 o 2 

The equations (I1.3P are simply the equations of the flow of H v \ M - see §2.4[ They are 
easily seen to imply (11.51) but some h corrections are built into the classical motion leading 
to the improvement in Theorem (TJ see also Fig{TJ 

As in previous works all of this hinges on the proximity of u(x, t) to the manifold of 
solitons, M. In [9] we followed 0] and used Weinstein's Lyapunov function [T3] . 

L(w) = f S(w + rf) - S(rj) , 
where £ is given by (11.81) . Here, in the notation similar to (11.91) . 

u(t)=g(t)-(r) + w(t)), 
for an optimally chosen g(t) = (a(t),v(t),j(t), (i{t)) - see Lemma [2731 

The use of L(w) seems essential for the all-time orbital stability of solitons. Up to times 
h~ x log(l//i) we found that it is easier to use its quadratic approxition 

(1.12) L = {£w,w), Cw = 8^(w) = -^d 2 x w-2r ] 2 w-r ] 2 w + ^w. 

Rather than use conservation of energy, Hy(u), we use the nonlinear equation for w(t) 
in estimating L Q (w) - see §|5j That involves solving a nonhomogeneous linear equation 
approximately using the spectral properties of C - see (15. 21) and Proposition 14.21 The 
fact that the solution of that equation is of size h 2 in H l gives the first indication of 
the improvement based on using the effective Hamiltonian. Ultimately, this makes the 
argument simpler than in the case of W G H' 1 (or the special case of W = 5q). 
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The paper is organized as follows. In £J2] we recall the Hamiltonian structure of the 
nonlinear flow of (11 .ip and describe the manifold of solitons. As in [9], its identification 
with the Lie group G = if 3 x IR+, where H3 is the Heisenberg group, provides useful 
notational shortcuts. In £J3] we describe the reparametrized evolution. The starting point 
there is an application of the implicit function theorem and a decomposition of the solution 
into symplectically orthogonal components. That method has a long tradition in soliton 
stability and we learned it from [I]. The analysis of the orthogonal component using an 
approximate solution to a linear equation and a bootstrap argument are presented in §5j 
This results in a somewhat more precise version of Theorem [1] - see Theorem |2j The 
ODE estimates needed for the exact evolution (II. 3p are given in §|6j Finally we show how 
Theorem [2] implies Theorem [TJ Except for basic material such as properties of Sobolev 
spaces or elementary symplectic geometry, and a reference to the proof of Proposition 14.11 
the paper is meant to be self contained. 

Acknowledgments. The work of the first author was supported in part by an NSF 
postdoctoral fellowship, and that of the second second author by the NSF grant DMS- 



2. The Hamiltonian structure and the manifold of solitons 

In this section we recall the well known facts about the Hamiltonian structure of the 
nonlinear Schrodinger equation. The manifold of solitons is given as an orbit of a semidirect 
product of the Heisenberg group and M + . 

2.1. Symplectic structure. In our work, we consider 



viewed as a real Hilbert space. The inner product and the symplectic form are given by 



Let H : V — > M be a function, a Hamiltonian. The associated Hamiltonian vector field is a 
map Eh '■ V — > TV, which means that for a particular point u G V, we have (Eh) u £ T U V. 
The vector field Eh is defined by the relation 



0654436. 



V = i? 1 (E, C) C L 2 (M,C) 



(2.1) 





(2.3) dH u (v) = (dH u ,v), (E H ) u = -dH u . 
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If we take V = H 1 ^, C) with the symplectic form (12. ip . and 

H(u) = J\\d x u\ 2 -\\u\" 

then we can compute 

d u H(v) = Re J {(l/2)d x ud x v - \u\ 2 uv) 
= Re J {-{l/2)d 2 x u-\u\ 2 u)v. 

Thus, in view of (Q and Q , 
The flow associated to this vector field (Hamiltonian flow) is 

(2.4) « = (S^)„ = | - i«r«) • 

For future reference we state two general lemmas of symplectic geometry. The simple 
proofs can be found in 0, §2]. 

Lemma 2.1. Suppose that g : V — > V is a diffeomorphism such that g*u = fi(g)u, where 
H(g) E C°°(V;R). Then for f e C°°(V,R), 

(2-5) (g-%E f (g(p)) = ^y^ r/ (p) , p E V . 

Suppose that / G C°°(V;M) and that df(p ) = 0. Then the Hessian of / at p , f"(po) '■ 
T p V i — > T*V, is well defined. We can identify T p V with T*V using the innner product, and 
define the Hamiltonian map F : T p V — > T p V by 

(2.6) F = -f"(p ) , (f(po)X, Y) = u(Y, FX) . 

% 

In this notation we have 

Lemma 2.2. Suppose that N C V is a finite dimensional symplectic submanifold of V, 
and f E C°°(V,R) satisfies 

~ f (p)eT p NcT p V, P eN. 

If at p £ N, df(po) = 0, then the Hamiltonian map defined by (12. 6p satisfies 

F(T P N)CT P N. 
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The Hamiltonian map, F, is simply the linearization of at a critical point of /. An 
example relevant to this paper is 

f(u)=S(u) d ^ H (u) + ^\\u\\ 2 , 

see (11. 8p . The soliton 77 is a critical point of £ and the Hessian of E is given by C in 
(j!.12p . The Hamiltonian map F = (l/i)C is the linearization of at 77. In other words, 
(l/i)(C — 1/2) is the linearization of (11. ip (with V = 0) at 77. The 1/2 term comes from 
||m|| 2 /4 in the definition of S. 

In Lemma [2.21 we can take iV to be the four dimensional manifold of solitons and p — rj. 
It then says that (l/i)£ preserves the symplectic orthogonality of w G T V V to T V M. 

2.2. Manifold of solitons as an orbit of a group. For g = (a, 7j,7,/i) G M 3 x M + we 
define the following map 

(2.7) H x 3u\ — >g-ueH\ (g ■ u){x) = e^e iv{x ~ a) fiu(ii(x - a)) . 

This action gives a group structure on M 3 x M+ and it is easy to check that this transfor- 
mation group is a semidirect product of the Heisenberg group H 3 and M + : 

G = H 3 kR +: n-{a,v,i) = {-,nv,i). 

Explicitly, the group law on G is given by 

(a, v, 7, fi) ■ (a', v', 7', //) = (a", v", 7", //) , 

where 

v" = v + v'/j, a" = a + — , 7" = 7 + Y + — , n" = /ip 1 

H /i 

The action of G is not symplectic but it is conformally symplectic in the sense that 

(2.8) g*u = fi(g)u, g = (h(g) , fi(g)) , fi(g) G R+ , 
as is easily seen from (12. ip . 

The Lie algebra of G, denoted by q, is generated by ei, e 2 , e 3 , e 4 , 
exp(tei) = (t, 0, 0, 1) , exp(te 2 ) = (0, t, 0, 1) , 
exp(te 3 ) = (0, 0, t, 1) , exp(te 4 ) = (0, 0, 0, e*) , 
and the bracket acts as follows: 

(2.9) [ei, e 4 ] = e u [e 2 , e 4 ] = -e 2 , [e 1: e 2 ] = -e 3 , [e 3 , •] = , 

so e 3 is in the center. The infinitesimal representation obtained from (12.71) is given by 

(2.10) e 1 = -d x , e 2 = ix, e 3 = % , e A = d x -x. 

It acts, for instance on S(M) C H 1 , and by X G Q we will denote a linear combination of 
the operators ej. 
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The proof of the following standard lemma can be again found in [91 §2]: 

Lemma 2.3. Suppose E 3 t \— > g(t) is a C 1 function and that u G <S(R). Then, in the 
notation of ( 12. 7p . 

j t g(t)-u = g(t)-(X(t)u), 

where X(t) G g is given by 

(2.11) X(t) = a(t)n(t) ei + 4|e 2 + (j(t) - a(t)v(t))e 3 + ^-e A , 
where g{t) = (a(t), v(t), -y(t), /i(t)). 

The manifold of solitons is an orbit of this group, G ■ rj, to which S#, defined in ( 12.21) . is 
tangent. In view of (12.41) that means that 

i {\ d l r i + l^l 2 ^) = x -v, 

for some X G 0. The simplest choice is given by taking X = Xi, A G M, so that t] solves a 
nonlinear elliptic equation 

~\ r l" - ^ 3 + A?7 = . 

This has a solution in H l if A = yU 2 /2 > and it then is r](x) = /isech(yux). We will fix 
/i = 1 so that 

7]{x) = sechx . 

Using Lemma [2.11 we can check that G ■ rj is the only orbit of G to which H# is tangent. 
We define the submanifold of solitons, M C Hi, as the orbit of r) under G, 

M = G ■ 7] C Hi 

and thus we have the identifications 

(2.12) M = G-rj~G/Z, T^M = g • rj ~ g . 
The quotient corresponds to the Z-action 

(a, v, 7, /i) i— > (a, 7 + 27r/c, /x) , k G Z 

2.3. Symplectic structure on the manifold of solitons. We first compute the sym- 
plectic form uj\m on T^M using the identification (12.121) : 



(u\ M ) v (ei, ej) = Im J (d ■ rj)(x)(e j ■ rf)(x) . 

Since 

r] 2 (x)dx = 2, / r](x)d x r](x) = , / d x r](x)xr](x)dx = — 1 
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we obtain from f)2.10p that 

(2.13) (^U/)r?(e 2 , ei) = 1 , (uj\ M ) v (e 3 , e 4 ) = 1 , 
and all the other (uj\M)r){ e i, e j)' s vanish. In other words, 

{u\m)t, = (dv A da + A dji) (0,0,0,1) = (d(uc/a + 7^))(o,o,o,i) • 
Using ( 12.121) we conclude that 

(2.14) w\m= fJ'dv Ada + vdfi A da + dj A d\i , 
see [H Lemma 2.3]. 

Now let / be a function defined on M, / = /(a, t>, 7, /i). The associated Hamiltonian 
vectorfield, Sy, is defined by 

S/) = d/ = / a da + + fudfj, + f^ . 

Using (I2.14p we obtain 

(2.15) ~ f = ^d a + (-h: - V Jl)d v + f 7 d, + (vh - A d, . 
The Hamilton flow is obtained by solving 

b 

The restriction of 
to M is given by computing by 

(2.16) f(a,v,j,fi) = H(g-ri) = — - — , g = (a, v, 7, . 

2 

The flow of ( 12.15ft for this / describes the evolution of a soliton. 

2.4. The Gross-Pitaevski Hamiltonian restricted to the manifold of solitons. We 

now consider the Gross-Pitaevski Hamiltonian, 

(2.17) ff ff («) = ^ J(\d x u\ 2 -\u\ 4 )dx + ± J V(x)\u\ 2 dx, 
and its restriction to M = G ■ 77: 

(2.18) H q \ M = f(a,v n ,fi) = ^--^ + ^V* (sech 2 (H)(a) • 



fa ./«.„. fv n 

, a = — , /i = / 7 , 7 = u 

H 11 /i /i 



H{u) = - f \d x u\ 2 -\ I 
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The flow of (H q )\ M can be read off from (12. 15^ : 
(2.19) 

fa Vf /i 2 /„ 

v = — = V' * (sech 2 (/j,«))(a) , a = — = v , fi = L = , 

fv 11 

7 = v— f fi — -v 2 + -/i 2 - fiV * (sech 2 (/j«))(a) + /j,V * (xsech 2 (a;) tanh.(x)\ x= ^.)(a) . 

This are the same equations as (jl.3p . The evolution of a and v is simply the Hamiltonian 
evolution of (v 2 + jj?V * sech 2 (yU«)(a))/2, /x = const. The more mysterious evolution of the 
phase 7 is now explained by (12.181) . 

We can also rewrite (12.191) as follows: 



so that 



and 



* = ~\ I v il + a ) v2{x) dx = ~ v ' {a) - v -^~ + o{h5) ' 

where we used that J rf = 2 and J x 2 rf = 7r 2 /6 in the Taylor expansions. Here for the 
purpose of presentation we assumed that W G C 5 but we will never use Taylor's formula 
with more than four terms, for which C 3 is only needed. 

3. Reparametrized evolution 

To see the effective dynamics described in §2.41 we write the solution of (11.11) as 

u{t) = g(t) ■ (77 + w{t)) , w(t) E H\R, C) , 

where w(t) satisfies 

u(w(t),Xri) = 0, Vleg. 
To see that this decomposition is possible, initially for small times, we apply the follow- 
ing consequence of the implicit function theorem and the nondegeneracy of u \m (see jH 
Proposition 5.1] for a more general statement): 

Lemma 3.1. For SdG/Z (where the topology on G/Z is given by the identification with 
R x R x S 1 x R + ) let 

Ut,,s = { u e Hi : inf \\u - g ■ rj\\ H i < 5} . 

gGY, 
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If $ < — ^o(S) then for any u G U^, there exists a unique g(u) G £ such that 

(3.1) ^(Vr 1 ■u-rj.X-r}) = VlGg. 
Moreover, the map u \— > g(u) is in C 1 (6 r s,5, £)• 

Proof. We define the transformation 

F : i7 1 (R,C)xG — ► g*, [F(u 7 h)](X)=u(h-u-ri,X-Ti), 

and want to solve F(u, h) = for /i = /i(u). By the implicit fuction theorem that follows 
for u near G ■ r\ if for any go & G the linear transformation 

4^0o • V, 9o) ■ T go G — ► g* , 

is invertible. Clearly we only need to check it for go = e, that is that dhF(rj, e) : — > 0* , 
is invertible. But as an element of 0* ®0*, dhF(rj, e) = (uj\m)ti: which is nondegenerate. □ 

From § §2.11 and 12.41 we recall that the equation for u (11. ip can be written as 

(3.2) d t u = E Hq (u) , H q (u) = i J (\d x u\ 2 - \u\ A )dx+ l - j V(x)\u\ 2 dx . 
Using Lemma [3. II we define 

(3.3) <?(*) = 9(u(t)), u^gitr'uit), w(t)^u- V , 

and we want to to derive an equation for w(t). 
Let 

a = a(a,jj) = - Vi — h a )r] 2 (x)dx / V' [ — h a ) — rfix) dx , 

( 34 ) 2 J \li J 2 J V/i / fi 



(3 = (3(a,fi) = ^- [v'(- + a)) V 2 (x)dx 



2 ix J \/J 

Note that dependence on a, ii makes a, (3 into time-dependent parameters. They are 
however independent of x. The proper motivation for the choice of a and (3 will come in 
Lemma 13.31 below. 

Set 

(3.5) X = (—a + v)iie x + ( - ~ - ^) e2 + ( — 7 + ^ v ~ 7^ 2 + \^ ~ a ) e 3 _ j^ e 4 , 

where e/s are given by 1 12. 10p . Let also 

(3.6) J\fw = 2\w\ 2 rj + r/w 2 + \w\ 2 w . 
We now have 
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Lemma 3.2. The equation for w , defined by f)3.3p . is 

V 

+ Xw + i - 

L v/i / 

where X is given by (j33| . C by (02]) . A/" 6y fM . 
Proof. We compute, by the chain rule 



(9(W = X77 + z — U ( — + a^j + a + j3x 



w 



ijji 2 Cw + ijj^M 



w 



where 



d t u = d t (g ■ (77 + w)) = g ■ Y(r] + w) + g ■ d t w 



r = a/iei + -e 2 + (7 — av)e 3 + -e 4 . 



From this, and the fact that the equation for u can be written dtu = we get 

d t w = —Yu + g~ l d t u = —Yu + g^^ngu 
We apply Lemma 12.11 to obtain 



d t w = -Yu + -E„* H u 



We compute 



(g*H)(u) = H(gu) = ^v 2 fi\\u\\ 2 L2 + ^fi 3 \\d x u\\ 2 L2 + ^vfi 2 lm \ 



-\^\H 4 L± + \^ I v(- + a))\u(x)\ 2 dx 



Therefore, 



, - 1/1 2 - 
^HU=-[-V^U 2 



c^u — ^ |u| u + fiVy — \- ajuj — vfi o x u 



Substituting and expanding the cubic term, 

d t w = — Y(r] + w) — -iv 2 (r] + w) + -i^ 2 d 2 {r] + w) + ifi 2 [r] 3 + 2r] 2 w + rfw 
+ 2|-u;| 2 ?7 + 7]w 2 + \w\ 2 w] — iV ^— + (77 + w) — vfid x (r] + w) 



Using that — \r\ + ^77" + 77 s = 0, we obtain 



fx \ 

dtw = — Yr] — —iv 2 r] + —i^ 2 r] — iVy — h ajr] — vfid x r] 

-ifi 2 w — ifi 2 Cw + iu 2 J\[w — iv( — + a)w — vfid x w 
2 \a J 



Yw iv 2 w , 

2 2 
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Now set 

X = —Y + v/iei — (3e2 
From this, we get the claimed equation. 
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1 2 1 2 

- r + r 



a 



e 3 



□ 



Let us make some remarks about the lemma. First, note that if X = 0, then 
v = —/3{i , 7 = av — -v 2 + -fi 2 — a , a — v , bfi — , 

which are exactly the equations of motion of the effective Hamiltonian - see §2.41 
Second, note that the term 



x 



V(- + a) +a + fix 



V 



projects symplectically to (used in the next lemma), in the sense that 

(3.7) P(i(-V(*/i2 + a) + a + f3*)r]) =0, 

where P : S'{R, C) — > q is defined by the condition that 

uj(u - P{u)rj, Yrj) = , VFGg 

see [H (3.7)] for an explicit expression of P. To see (13.71) we use the following simple 

Lemma 3.3. Let P be the symplectic projection defined above and f G «S'(R;R) (that is f 
is real valued tempered distribution) . Then, 

P(if{x)rj) = aez ■ rj + [3e2 ■ r\ = iar] + ij3xrj , 

(3 8) 1 f 1 f 1 f 

a= 2 f^rf^ ~ 2 / f'( x ) x V 2 (x)dx , P=2 f'{x)r] 2 {x)dx. 

Proof. This follows from a straightforward calculation based on (12.131) . If f(x) = V(x/fi + 
a), we obtain (13.71) with a and /3 given by (13. 8p . □ 

Finally, note the following Taylor expansions: 



■[ - + a) = V(a) + V'(a)- + V'\a)^ 
a = V(a) - + 0(h*) , 



0(h 3 







and thus 
(3.9) 



x 



-V[- + a)+a + (3x 



24/x 2 
+ 0(h 3 ) , 

V"(a)7r 2 V"t 



24/i 2 



2/i 2 



-x 



0(h 3 ) = 0{h 2 ) 
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where all the errors are polynomially bounded in x. If we assumed that W G C 4 then he 
expansions for a would be valid with an error 0(h 4 ). 

Lemma 3.4. Let w be given by ( 13.31) . with g obtained from Lemma \3. 11 and X by ( 13.51) . 

Suppose g < fi < 1. T/ien 

|X| < c{h 2 \\w\\ H i + ||w||^i + \\w\\%i) 

Proof. We use the symplectic orthogonality of Yrj, Y G and u>. 
Since Piwt = <9iPw; = 0, Lemma [3T21 gives 

X = P (i(V(m/fi + a )-a- P*)r)) + P(i(V(*/fi + a) - a - /3m)w) - P(Xw) 
+ n 2 P(iMw) + ii 2 P(iCw) . 
We recall from 0, Lemma 3.3] the following straightforward estimates: 

(3.10) \\P(Yw)\\<C\Y\\\w\\ L 2, \\P(iMu)\\<C\\w\\l 2 (l • 

We already observed in (13.71) that the first term on the right hand side vanishes^. From 
(13. 9p we see that the second term is C(/i 2 )||w||l2. The third and fourth term are estimated 
using (13.101) by 

C(\X\\\w\\ L 2 + \\w\\ 2 Hl + \\w\\ 3 Hl ). 

The last term vanishes: the linear operator L is the Hessian of S, given in (11.81) . at the 
critical point rj. The fact that is tangent to M and Lemma [2^2] (or a direct computation) 
show that 

P(iCw) = , 

Summarizing, 

\X\ < C||«;|| L 2|X| + C(h 2 \\w\\ H i + \\w\\ 2 H i + \\w\\ 3 H i) . 
The smallness of ||w||l2 concludes the proof. □ 

We conclude this section which two lemmas which effectively eliminate [i from the coef- 
ficients of X. 

Lemma 3.5. Suppose that w G PT 1 (1R, C) and that u(w, X ■ rj) = 0, for every leg. Then 

(3.11) H|| a = 2(1 ii. 

Proof. We first compute 

\\n + w\\ 2 L 2 = Wg^uW 2 ^ = \\u\\jp/n{g) = 2//i(g) , 

where we used the conservation of the L 2 norm. By the symplectic orthogonality assump- 
tion statement of the lemma Ke(w,r]) = and hence 

||t7 + HIz,2 =2 + \\w\\ 2 L 2, 
^This is the significant bonus of using the effective Hamiltonian. 
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□ 



From this we immediately deduce the following 
Lemma 3.6. Suppose l-/j< 1 and < h < 1. Let X = X\ fJL= i where X is given by 

(3.12) X d = (—a + v)e 1 + ^ — v — /^e 2 + - 7 + - ^v 2 + - - «je 3 - /ie 4 , 



where a and (3 are given by (13.41) (and depend on The the conclusions of Lemma 3Jj_ 
hold for X : 

\X Q \ < c(h 2 \\w\\ H i + \\w\\ 2 H i + \\w\\ 3 H i ) . 



4. Spectral estimates 

In this section we will recall the now standard estimates on the operator C which arises 
as Hessian of £ at 77: 

Cw = —-d 2 w — 2r] 2 w — r] 2 w + —w , 



or 



Cw 



~L + " 




Re w 


L_ 




Imto 



■^ = - 1 -d 2 x -(2±i) v 2 + 1 - 



In our special case we can be more precise than in the general case (see [H], and also [U 
Appendix D]). The self-adjoint operators L± belong to the class of Schrodinger operators 
with Poschl- Teller potentials and their spectra can be explicitly computed: 

<r(L_) = {0} U [1/2, 00) , <r(L+) = {0, -3/2} U [1/2, 00) . 

The eigenfuctions can computed by the same method but a straightforward verification is 
sufficient to see that 

3 , 



L„r/ = 0, L+(d x r}) = 0, L + ( V 2 



From §P, §4] we recall the following 

Proposition 4.1. Let w G -fT 1 (!R, C) and suppose that for any I 6 5, u(w,X • 77) = 0. 
Then, 

9 



(4.1) 



(£w, w) > 2p ° \\w\\ 2 m ~ 0.0555111011^! , p 



7 + 2p "^"" i r " 2(12 + vr 2 )- 

The next proposition will be useful in solving a linear equation for an approximation of 



w. 
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Proposition 4.2. The equation 

(4-2) L+f = + ^) ^) 

has a unique solution in L 2 (M). In addition, 

e {1 - e)lxl f {k \x) G L°°(M) , Ve>0, fceN, 

and 

(4.3) u(f,X v )=0, VlGg. 

Proof. Since is an isolated point of the spectrum of L + and L + is self-adjoint, L + has a 
bounded inverse on 

(kerL + ) ± = (span{9 x ?7})- L . 

Hence, 

V \ L- 1 



(ker^)- 1 9 + x 2 j r) i— ^ / G (kerL H 

By elliptic regularity, / is smooth. Moreover, since L + commutes with g(x) \-* g(—x) 
and + x 2 )rj(x) is even, we conclude that / is even. 

Next we argue that / has exponential decay. For that we conjugate the equation by e ax : 

e ax L +e - ax = L + + ad x - a 2 /2 , 

and apply both sides to e ax f(x): 

{^ + * 2 )v=(-dl + 2ad x - £ + i)e"/ + 3 V 2 e° x f . 
Taking the Fourier transform this means that we have 

which makes sense if \a\ < 1 as then the right hand side is in L 2 . Since the multiplier 
on the left-hand side is bounded away from 0, we can invert the expression to obtain that 
e ax f G H 2 . From this, we deduce by Sobolev embedding that e ax f G L°°. By applying 
derivatives to (14. 2[) and then repeating this argument, we in fact obtain that all derivatives 
are pointwise exponentially localized in space. 

Finally, we prove the orthogonality condition (I4.3p . that is, u>(f, Cj -rj) = where e/s are 
given by f)2.10p . Since / is real, we clearly have 

«;(/,»/)="(/, M') = 0. 



SOLITON INTERACTION WITH SLOWLY VARYING POTENTIALS 



17 



Since / is even, we also have cu(f,ixr]) = 0. It remains to show that u(f,ir]) = 0, that is 
that j ft] = 0. Note that L + [(xf]) f } = 77 (by direct computation or from the structure of 
the generalized kernel of iC). Hence 

Integration by parts and j rj 2 = 2, j x 2 rf = n 2 /6, show that this is 0. □ 

5. Proof of the main estimate 



In the arguments that follow, we will assume that g(u) satisfying (13.11) is always defined 
on any time interval under consideration, and thus w given by (13.31) is also always defined. 
This is, however, not known a priori and Lemma 13.11 must be considered as part of the 
bootstrap argument, together with the lemmas that follow. However, since this aspect of 
the argument is standard in papers on this subject, we will not make mention of it. 

Recall from Lemma [3.21 that the equation for w is 



(5.1) 



d t w = Xrj + i —V[ — \- a) + a + j3x 



x 



+ Xw 



V[- + a 



a + (3x 



w — ifi 2 Cw + ifi 2 Afw 



where a and (3 are time dependent parameters given by (13.41) . Note that by Taylor's 
theorem the forcing term in this equation has second-order expansion 



— 1 



x 



VI - + a) + a + (3x 



V"(a) (tx 2 



+ x 2 )r] + h 3 (x) 3 r)f 2 



2/i 2 U2 

where, provided! 1/2 < /x < 1, we have |/ 2 (a;, i)| < c. 

We first consider the forced linear evolution obtained from (15. ip by discarding all terms 
we expect will be of order h 3 or higher: 



(5.2) 



d t w = —i[i 2 £w 



Y"(a) /7T 



+ x }rj 



2/i 2 V12 

We now introduce a natural approximate solution to this forced linear equation. Let 

V"(a) 



(5.3) 



w 



where / is given in Proposition 14.21 Then w satisfies (15. 2p to second-order, i.e. 
(5.4) - ' V " M ^ 



d t w = -ifCw - i ^(^ + x 2 )ri + hHS 



■i-This will follow easily from the bootstrap assumption. 
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where 



def 



1 



V"\a)a 2V"(a)fi 



2/i 4 /i 5 
Lemma 5.1. Suppose there is a constant c\ such that 
(5-5) . ,Hl < C!h 3/2 

Then provided \t 2 — t±\ < h~ 2 , we have 

sup \9(t)\ < c, 

t!<t<t 2 

where c is a constant depending only upon c\, 

IIW^IU- fork = 0,1,2,3 and |v(t x )|. 

Proof. By (15. 5p . Lemma [3.61 and Taylor expansions (see (13.91) ). we have 

(5.6) |u + y'(a)| + + | -a + v\ < ch 3 

By integrating the first bound, we obtain the following rough estimate 

(5.7) sup |u(*)| < + c/illW'Hioo^a — *i) + c/i 3 (* 2 — *i) ■ 

t!<t<t 2 

We have 

\vv + V'(a)a\ <\v + V'(a)\\v\ + \V'(a)\\a - v\ 
and thus, by (15. 6p and (15.71) . 

sup \vv + V'(a)a\ < c/i 3 |t>(ti)| + ch 4 \\W'\\L°°{t 2 - h) + ch 6 {t 2 - h) 

tl<t<*2 

Integrating this bound, we obtain a near conservation of classical energy, 

(j+na))-(^f + v {a{k) ))\ 

< c/i>(*i)|(t 2 - ti) + ch A \\W'\\ L °°(t 2 - ti)(t 2 - ti) + ch 6 (t 2 - h) 2 
from which we obtain sup tl<t<t2 \v(t)\ < c. This and (I5.6P imply that 

sup \a(t)\ < c, 

ti<t<t 2 

and thus \6(t)\ < c. □ 

This approximate solution w provides heuristic evidence that w should be of order h 2 , 
but it will also play a key role in our rigorous argument establishing this fact. 

Lemma 5.2 (Lyapunov energy estimate). Suppose that, for some constant c±, 

(5-8) \\w\\ L ~ , M <C!/i 3/2 

Then, provided 

(5.9) |* 2 -*i|<p 
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we have, 



(5.10) |HU~ , ,w < 4 v /^||u;(ti)||i f i + ch 2 

l*ls*2J 

The constants c in (15. 9ft and fl5.1Up depend upon c\ } ||W( fe )|| L oc for k = 0,1,2,3, and 

although are independent of . The constant C2 = (7 + 2po)/(2po) ~ 18.02, with 

po the absolute constant appearing in (14.10 . 

We will postpone the proof of Lemma 15.21 to the end of the section. In the next theorem 
we iterate the above bound, and close the bootstrap argument. 

Theorem 2. Let < 5 < \. Let Wq = w(0), and suppose that there is a constant c\ such 
that 

||w ||ffi < Ci/l2 +3<5 . 

Then, provided 

c(l + 5|log/i|) 

(5.11) |t| < ^ andO<h<e, 
we have 

(5.12) HHIi^ < ^h- 2S \\w \\ m + ch 2 ^ 

The constant c in (15. lip and (15.121) and the constant e in (15.111) depend upon c\, WW^Wlcc 
for k = 0,1,2,3, and \v(0)\, although are independent of 5 and \\wqWh 1 - The constant 

c 2 = f (7 + 2p )/(2po) ~ 18.02, with p the absolute constant appearing in ( 14. II) . 

Proof. We apply Lemma [5.21 k times on successive intervals each of size c/h (where c is as 
given in Lemma 15. 21) to obtain the bound 

\\w\\L~ ch/h] m < (4v^ 2 ) fe |KI|^ + ^(V^'j ch 2 

This is only valid provided that the hypothesis of Lemma 15.21 is satisfied over the whole 
collection of time intervals: 

(V^kii* + (j2^v~c 2 y^j ch 2 < Cl h^ 2 , 



25\\ogh\ 



By taking 

k = 

we obtain that 

(k-l 
J2^v~c 2 y ) ch 2 < (A^h^wwoW^ + ch 2 ^ 
j=0 
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Thus, it suffices to ensure that 

(4v^ 2 )/r 25 |KllHi + ch 2 ^ < Cl 2 , 

and this is accomplished provided h < e, for a suitable e with the dependence as stated in 
the proposition. We note that the constant c from Lemma 15.21 did not change from one 
iteration to the next, since \v remains uniformly bounded by Lemma |5"TT| which applies 
on a time interval of size h~ 2 . □ 



We now conclude the section with the proof of Lemma | 
Proof of Lemma \5. 6 A For the proof, we shall assume, in place of (15.81) . the bound 

IMUj? t & < ci/i 2(1 - 5,) 

1*1 .*2 J 

We will conclude at the end that it suffices to take 8' — ~. Let 

def 

W\ = W — W , 

where w is an approximate solution to (15.21) given by (15.31) . 
From ( 15.11) and ( 15 .4ft . we derive the equation satisfied by w\\ 



d t wi 



ifi Cw\ + Xrj + i — V [ — h a I + a + fix + 



x 



V"(a) (tx 2 



2{j? V12 



— + x' 



7] - h 3 6f 



+ Xw + i 



X 



-V[- + a)+a J r(3x 



w + in 2 Nw . 



By grouping forcing terms of size 0(h 3 ), we rewrite the above as 

d t wi = —ifi 2 Cw\ + h 3 fi + Xw + i — v(^— + a^j + a + /3x w + ifi 2 J\fw 

where t)\\ip. < c (uniformly in t). Note that here we have applied Lemma 15.11 to 

conclude that \9(t)\ < c. 

We recall that £ is self-adjoint with respect to 

(u, v) = Re / uv , 



and hence, writing Xw = Xw\ + Xw, we compute 
-d t {Cw u w x ) = (£w u d t wi) 



/j, 2 (£w±,i£wi) + (£wi, h 3 fx) + (£wi, Xwi) + (£wi, Xw) 



+ (Cwi, i 



x 



V[- + a) +a + fix 



Wi) 



+ (£wi, i — v(^— + aj + a + (3x w) + fi 2 (£wi,iJ\fw) 



I + II + III + IV + V + VI + VII 
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Now we analyse these terms one-by-one. First, 

(5.13) 1 = 0. 

For II, we use integration by parts for the d 2 term to move one d x onto fx and then 
apply the Cauchy-Schwarz inequality; for the other terms we use a direct application of the 
Cauchy-Schwarz inequality, and together these give 

(5.14) | II | < ^\\wx\\ H 4fx\\ Hl < ch 3 \\wx\\ H i 
The next term, III, requires more care: 

III = (£wx,Xwx) 

= _ dl w i - 4?7 2 u>i - 2i] 2 wx), (-axd x wx + a 2 ixwx + a 3 iwx + a±d x (xwx))) 

where cij, the components of X, are time dependent but space independent. By the boot- 
strap assumption (15.81) and Lemma 13.41 \aj\ < ch 4 ~ iS . To proceed, we do some further 
calculations: 

(wx,Xwx) = a 4 (wx, (wx + xd x Wx)) = — a 4 ||^i Ilia 

(d^wx, Xwi) = a 2 (dlwx,ixwx) + a^dlwi.wx + xd x wx) 

3 

= a 2 (d x wx,iwx) - ^a 4 \\d x Wx\\ 2 L 2 

and thus the above two terms are bounded by c/i 4_4<5 ' ||u>i||^i. For the terms involving 77, 
we use the Cauchy-Schwarz inequality and the fact that rj £ S: 

\(r] 2 Wx,Xwx)\ + \(r] 2 Wx,Xwx)\ < (max |aj|)||wi||^i 
Altogether then, we have 

(5.15) |III| < ch 4 - u '\\wxf m 
Now we move on to IV: 

IV = (Cwx,Xw) 

= -((wx — d 2 w\ — Arfwi — 2r] 2 wx), {—axd x w + a 2 ixw + a 3 iw + a i d x {xw))) . 

For this term, we are forced to directly estimate by the Cauchy-Schwarz inequality (al- 
though for the d 2 Wx term, we integrate by parts one d x factor). The bound that we get 
is 

|IV| < ch 4 - 45 '\\wx\\m(\\(xM\ L z + \\(x)d x w\\ L 2 + || (a;>«9 x 2 w;|| L2 ) 

From the definition (15.31) of w and Proposition 14.21 we obtain that all the norms involving 
w are bounded by h 2 . Thus, 

(5.16) |IV| < ch 6 ~ 4S ' \\wx\\ m ■ 
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Next, we move on to V: 



V = (Cwi, i — V ^— + a J + a + (3x Wi) 



— {{w\ — d 2 Wi — 4r] 2 wi — 2r] 2 wi),i — v(— + a) + a + fix Wi) 



x 



-V[ — + a) + a + fix 



Wi) 



< ch 



To estimate the first term, we integrate by parts and use that 

/i V/i / 

(note that other estimates are available, like cxh 2 , but we do not want an x coefficient 
here). For the second term, we use f !3 . 9H . the Cauchy-Schwarz inequality and the rapid 
decay of rj 2 : 



-V[- + a) +a + fix 



<ch 2 



This gives the bound! 
(5.17) \V\ < ch\\ Wl \\ 2 Hl . 

Now we move on to the next term, VI. 

VL = {Cw u i -V^- + aj +a + fix w) 

In the d 2 term of C, we integrate by parts one d x , and then estimate by the Cauchy-Schwarz 
inequality. All other terms, are estimated by a direct application of the Cauchy-Schwarz 
inequality. The bound obtained is 



I VI I < llwillffi 



(d x ) 










4- a + (3x 


w 









L' 2 



< ch 2 \\wi\\ H i(\\(x) 2 w\\ L 2 + \\(x) 2 d x w\\ L 2) 
(5.18) < ch 4 \\wi \\ H i 

by the localization of w. For the last term, VII, we use integration by parts once for 
the d 2 term, and then apply the Cauchy-Schwarz inequality to all terms. Since we are in 
one-dimension, we have the embedding ||iu||z,°o < c||iu||#i. 

-fi 2 (£wi, iAfw) 



(5.19) 



by the bootstrap assumption. 



VII 



|VII| < + \H\m) < c/i^'IKIIh 1 



^It is unlikely that we can do better than h as a coefficient here, and thus this seems to be what limits 
us ultimately to time 1/h. 
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This completes the step- by-step estimation process, and the bound we get from (I5.13p . 

(EH, (JEH, ALU, dEIZD- dHHD, dEH is 

\d t (Cw 1 ,w 1 )\ < c(h^ 5 ' + h 3 )\\ Wl \\ m + c(h 4 ^ s ' + h)\\ Wl \\ 2 Hl 
We see that it suffices to take 5' = \. Integrating in time, we get 

(5.20) (Cwi(t),wi(t)) < (£w 1 (t 1 ),w 1 (t 1 )) + c(t-t 1 )/i 3 ||u ;i || i cx ) H x 

L* 1 -*2J 

+ c{t - ti)h\\wi\\l°° t , H i 

L*1>*2J 

By (13. ip and (14.31) . Wi(t) satisfies the hypothesis of Proposition 14. 1[ and we have 

(5.21) -HWIIiri < (£wi(t),wi(t)) 

C 2 

By direct estimation, we have the upper bound 

|(£wi(*i),wi(*i))| < 411^(^)11^1 
Combining this with (15.201) we get the bound 

IN0)||tfi < 4c 2 ||u7i(ti)||^i + c(t - ti)/i 3 |K|| L ~ , ,m 
+ c{t -ti)/i||wi||^ Hl 

[*l.*2j 

From this, we infer from the monotonicity of the right side that 

< 4c 2 ||wi(ti)||^i + c(t 2 - ^i)^ 3 ||wi||l- , ,hi 

[i 1 ,i 2 ] 1*1. *2J 

[*1.*2] 

Requiring that t 2 — ti < c/h implies 

Ikilli- < 8c 2 ||wi(ti)||^i + c/i 2 ||wi|| Lf » m 

[ti,*al 1*1.421 

< 16c 2 ||wi(ti)||^i +ch A . 
Since w = W\ + w and ||w||//i < ch 2 , 



(5.22) ||w|| L jx, „i < 4^11^(^)11^1 + ch 2 , 

which is the claimed estimate. □ 



6. ODE ANALYSIS 

To pass from an approximate equations for the parameters of the soliton, (a, v, 7, /z) given 
in Lemma 13. 6} to the ODEs (II. 3p we need some elementary estimates which we present in 
this section. They are similar to those in [HJ §7]. 



(6.1) 
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Lemma 6.1. Suppose that < h <C 1, and a = a(t), v = v(t), €\ = e\(t), e 2 = ^{t) are 
C 1 real-valued functions . Suppose f : R — > IR is a C 2 mapping such that \ f\ and \ f'\ are 
uniformly bounded. Suppose that on [0, T], 

a — v + ei a(0) = a 

v = hf(ha) + e 2 ' v(0) = v 

Let a = a(t) and v = v(t) be the C 1 real-valued functions satisfying the exact equations 

a = v a(0) = a 

v = hf(hd) ' v(0) = v 

with the same initial data. Suppose that on [0, T], we have \ej\ < h 4 ~ s for j = 1,2. Then 
provided T < 5h~ x log(l//i) ; we have on [0,T] the estimates 

\a-a\< h 2 ' 25 log(l//i), \v - v\ < h 3 ~ 25 log(l//i) 

Before proceeding to the proof, we recall some basic tools. 

Gronwall estimate. Suppose b = b(t) and w = w(t) are C l real- valued functions, h is a 
constant, and (b, w) satisfy the differential inequality: 

\b\ < \w\ b(0) = b 

\w\<h 2 \b\ : w(0)=w 
Let x{t) = hb(t/h), y{t) = w(t/h). Then 

x(0) = x = hb 
2/(0) = y = w 

Let z(t) = x 2 + y 2 . Then \z\ = \2xx + 2yy\ < 2\x\ \y\ + 2\x\ \y\ < 2(x 2 + y 2 ) = 2z, and hence 
z{t) < z{0)e 2t . Thus 

\x(t)\ < v2 max(|xo|, \yo\) exp(t) 
\y(t)\ < V2m&x(\x \,\y \)exp(t) 
Converting from (x,y) back to (b,w), we obtain the Gronwall estimate 

\b(t)\ < v / 2max(/i|6 |, \w \)^$- 

\w(t)\ < V2max(h\b \,\w \) exp(ht) 

Duhamel's formula. For a two- vector function X(t) : IR — > M 2 , a two- vector X G M 2 , and a 
2x2 matrix function A(t) : R — ► (2 x 2 matrices), let ^(t) = 5*(t, t')X denote the solution 
to the ODE system X(t) = A(t)X(t) with X(t') = X : 

j t S(t, t')X = A(t)S(t, t')X , S(t', t')X = X . 



(6.2) 



r \x\ 


< 


\y\ 


I \y\ 


< 


\x\ 
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— > IR 2 , the solution to the inhomogeneous 



Then, for a given two- vector function f(t) : 
ODE system 



(6.4) 



X(t) = A(t)X(t) + F(t) 



with initial condition X(0) = is given by Duhamel's formula 



(6.5) 



X(t) 



S(t } t')F(t')dt' 



Proof of Lemma \6. 1\ Let a = a — a and v = v — v ; these perturbative functions satisfy 

a = v + €\ a(0) = 

t = h 2 ga + e 2 ' ^(0) = 



where g — g(t) is given by 

f(ha) - f(ha) 
h(a — a) 
f'(ha) 

which is C 1 (in particular, uniformly bounded). Set 



if a ^ a 
if a = a 



A(t) 



1 

h 2 g(t) 



F(t) 



e 2 (t) 



X(t) 



~a(t) 

m 



in (16.41) . and appeal to Duhamel's formula (16.51) to obtain 



(6.6) 



S(t,t>) 



e2(t') 



dt' 



Apply the Gronwall estimate (16.31) with 

S(t + t\t') 



b(t) 
w(t) 







'b ' 




'exit')' 




5 


_w _ 







to conclude that 

S(t,t') 



/i" 1 exp(/i(t - t')) 
exp(h(t-t')) 



max(/i|e 1 (t , )|,|e 2 (OI) 



Feed this into (16.61) to obtain that on [0,T] 

\a(t)\ < v^T 6 ^^ sup max(/z|ei(s)|, |e 2 (s)|) 

h 0<s<T 

\v(t)\ < V2Texp(hT) sup max(/i|e 1 (s)|, |e 2 (s)|) 

0<s<T 

Taking T < 5h~ l log(l//i), we obtain the claimed bounds. 



□ 
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7. Proof of Theorem Q] 

We can now put all the components of the proof together. Lemma 13.61 and Theorem [2] 
show that on the time interval < t < c5\og(l/h)/h we have (II. 2p with the parameters 
satisfying 

d = v + 0(h^- s) ) , v = -sech 2 *V' / (a)/2 + 0(/i 4(1 - <5) ), fi = 0(h^- &) ) , 
7 = 1/2 + v 2 /2 - sech 2 * V(a) + (x sech 2 xtanha;) * V{a) + C(/i 4{1 ~ 5) ) . 

Lemma 16.11 can be applied to replace a and v with solutions of ( 11 .3D and the direct inte- 
gration of the error terms shows that the same is true for \x and 7. In particular we can 
drop n altogether. □ 

We conclude the paper with some remarks. The proof above and Theorem [2] show that 
the conclusions of Theorem [1] remain unchanged if instead of taking e lx1 ' sech(x — Xq) as 
initial condition, we took 

e ixv °sech(x - x ) + r(x) , ||r||#i < Ch 2 ~ s . 

We could go down to < Ch 3 ^ 2+3S at the expense of complicating the final statement 

to (I5.12p . A more general condition on the initial value would make the bootstrap argument 
in §5] so unwieldy that we opted out of pursuing that technical issue. 

In higher dimensions similar methods are clearly applicable for weaker nonlinearities and 
under further spectral assumptions - see [4] for examples. At this early stage we restrict 
ourselves to the physically relevant cubic nonlinearity which at the moment is tracktable 
only in dimension one. 
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